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Abstract: The supersymmetric SO(10) GUT with t-b-r Yukawa coupling unification 
has problems with correct electroweak symmetry breaking, experimental constraints 
(especially b —>■ 37) and neutralino abundance, if the scalar masses are universal at 
the GUT scale. We point out that non-universality of the scalar masses at the GUT 
scale generated both by (1) renormalization group running from the Planck scale to 
the GUT scale and (2) D-term contribution induced by the reduction of the rank 
of the gauge group, has a desirable pattern to make the model phemenologically 
viable (in fact the only one which is consistent with experimental and cosmological 
constraints). At the same time the top quark mass has to be either close to its quasi 
IR-fixed point value or below ~170 GeV. We also briefly discuss the spectrum of 
superpartners which is then obtained. 
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Grand Unification (GUT) has been regarded as a serious candidate of physics 
beyond the weak-scale to explain various puzzling features of the standard model. 
While the three gauge coupling constants fail to unify in its simplest version, they 
meet at a scale Mqut — 2 x 10 16 GeV in its supersymmetric (SUSY) extension. 
SUSY-SO(IO) GUT offers further exciting possibility that all three Yukawa coupling 
constants of top, bottom quarks and tau lepton may also unify at the same scale 
where the gauge coupling constants unify. This is possible because the supersymmet- 
ric standard model contains two Higgs doublets, and large values of tan/3 = v%/v\ 
(the ratio of the vacuum expectation values for the doublet H 2 and Hi which couple 
to the up and down quarks, respectively) lead to a proper bottom top quark mass 
hierachy, with approximately equal b and t Yukawa couplings [|TJ . The consequence 
of such an exact unification of couplings is that the top quark mass, nit, and tan/? 
are determined, once the bottom quark mass, nib, the tau lepton mass, m T , and the 
strong gauge coupling, a^, are fixed [|], ||, [|J. 

In this context, an interesting question is the issue of the compatilibity of this 
exact Yukawa coupling unification with the possibility of breaking the electroweak 
gauge symmetry through radiative effects. This question has been investigated in 
a number of papers in the minimal SUSY-SO(IO) models with universal IJ, ||, 
|| and non- universal [7]]-[]10|| soft supersymmetry breaking parameters at the GUT 
scale. Moreover, it has been recently observed that for these large values of tan (5, 
potentially large corrections to m;, may be induced through the supersymmetry 
breaking sector of the theory [||, ||. Altogether, the requirement of a physically 
acceptable value for the m& and of the consistency with the recent CLEO result for 
the decay b —>■ s'j and with the condition Qh 2 < 1 for the relic abundance of the 
LSP strongly constrain the minimal SUSY-SO(IO) with radiative electroweak sym- 
metry breaking. In recent papers |[L0|| , [|IT| it has been shown that those constraints 
rule out the model with universal soft supersymmetry breaking terms at the GUT 
scale and select certain class of non-universal boundary conditions which lead to 
radiative breaking with M 2 S> /i, i.e. with higgsino-like lightest neutralino, as the 
only acceptable scenario for the minimal SUSY-SO(IO). 

From the theoretical point of view, non-universal SUSY breaking terms at the 
GUT scale appear at present as a realistic possibility. In GUT models, even with 
universal boundary conditions at the Planck scale, the renormalization group (RG) 
running to the GUT scale generically leads to some non-universality of the scalar 
masses at that scale |§, fll2| . In addition, in models like SO(10), the reduction of 
the rank of the gauge group by one at Mqut-, together with non- universal scalar 
masses, generates additional non-universal contributions given by the D-term of the 
broken U(l) 0,0. 

In this paper, motivated by the phenomenological analysis of ref. |10|, we point 
out that the combination of both types of effects in the minimal SUSY-SO(IO) nat- 
urally gives the physically desirable non-universal boundary conditions at the GUT 
scale. This is a non-trivial prediction of the minimal model which depends in a 
crucial way on the presence of the D-term contribution (but not on its actual value) 
and on the value of the top quark Yukawa coupling, h t . Two branches of correct 
solutions are obtained: one for h t very close to its quasi- infrared fixed point and 
the other one for lower values of h t . With the present uncertainty on the top quark 
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mass, both solutions may be of phenomenological interest. 
The Higgs potential 

V = m\H\Hi + m\H\Hi — rn\ {H\iTiiii + h.c.j + quartic terms 



(m 2 = mjj. + /i 2 , m| = Bfi where \x is the supersymmetric Higgs mixing parameter 
and B is the corresponding soft term) has for large tan f3 values two characteristic 
features. It follows from the minimization conditions that 

o Ml 

m\ ~ —± (2) 

and 

9 Ml 

with 

Mj| ~ ml + m\ > 0. (4) 

Equations (0) and (|3|) are the two main constraints on the parameters of the scalar 
potential, which are characteristic for large tan f3 solutions. Combining @) and (gj) 

we get 

m\-m\> Mf . (5) 

Let us first discuss the dependence of the low-energy parameters on the GUT- 
scale boundary values. This discussion clearly shows the need for a special type 
of non-universality of scalar masses at the GUT scale. After that, we demonstrate 
such a special type of non-universality can be naturally obtained in SO (10) GUT 
even with a universal boundary condition at the Planck scale. 

The parameters of the low energy potential are given in terms of their boundary 
values at large scale and the RG running. Here we consider the minimal SUSY- 
SO(10) model with matter fields in 16 dimensional representations and the two 
Higgs doublets in one representation of dimension 10. At the GUT scale, after RG 
running from the reduced Planck scale Mp\ ~ 2.4 ■ 10 18 GeV to the GUT scale with 
the SO(10) RG equations, all scalar masses depend on three parameters 
and D (D-term): 

m Hi,H 2 = { _o \ D i (6) 




™ 2 q,v,d = ™ 2 m+{ 1 }D. (7) 

The values of the masses at the electroweak scale are obtained by solving the RG 
equations of the minimal supersymmetric standard model with the initial condi- 
tons at the GUT scale given by eqs.(|) and (0). For Y t = Y b = Y = h 2 /Air the 
approximate solutions to the 1-loop RG equations read [[|: 

m H lt H 2 = (l - jVzj mlo ~ ^yzm 2 16 - c M M 2 1/2 + j _^ | D + . . . , (8) 



™ 2 q,u,d = ( 1 " \vz) rnle ~ + <W> D mI /2 + { 1 \ D + ... . (9) 
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Mi/2 is the gaugino mass at M GUT , cm ~ C(2), c M ,U ' D ~ 0(A) and for convenience 
of analitical solutions we have introduced the parameter 

», - ™ (10) 

where Y (Mz) in the well known solution to the MSSM renormalization group equa- 



tions for the Yukawa coupling in the limit of large tan/5 [14 



y{ )= E lISSM (M 2 )Y(Mg VT ) 
y ' 1 + 7F mssm (M z )Y(Mcvt) y 

Here, Yfz is the auxiliary parameter given by the value of Y (Mz) corresponding to 
the Landau pole in Y at M GUT 



Emssm (Mz 
Y(Mgut)-^oo~ 7F MSS m(Mz 



Y fz = „,„ lm i Y \ M z) = ^ 77TT ' ( 12 ) 



and Emssm and Fmssm are functions of the gauge couplings. The dots in eqs.(|[) 
and (H) stand for terms which depend on soft parameter A at the GUT scale. In 
this approximation the condition (|5|) gives 

ml - m 2 2 = 4D > Mf . (13) 

However, here we have neglected small differences in the running of the two Higgs 
masses which follow from the different hypercharges of the right top and bottom 
squarks, from the difference in the running of the bottom and top Yukawa couplings 
(equal at the GUT scale) and from the effects due to the r lepton Yukawa. After 
inclusion of those effects we get 

ml - ml = aM 2 1/2 + cm\ + AD , (14) 

where m is an average scalar mass at the GUT scale (actually ml = (mf + 2mf 6 )/3) 
and the numerical values of the coefficients are a ~ |c| ~ 0(0.1) with c < 0. 
The small effects neglected in eq.([U|) but included in eq. (|T4]) are resposible for 
radiative breaking in the case of universal boundary conditions at the GUT scale 
(D — , mi 6 = mio = m ) 0. Then the large tan/3 solutions must be driven by 
large values of Mi/ 2 : 

M x /2 > ^ , Mi/a > m , (15) 

V a 



and as discussed in ref.jll], this scenario is strongly disfavoured for several reasons. 



It is clear from eq.flT4D that in the framework of SUSY-SO(IO) , with D = dm, 



qualitatively new solutions become possible if c + Ad > 0, with Mi/ 2 ~ and 

- ^ Mz n R \ 

mo> 7^- (16) 
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Thus, with positive D, contrary to the universal case, radiative electroweak breaking 
can be driven by soft scalar masses and this pattern does not depend on the actual 
value of the D term as well as on the values of mie and m w . However, as we shall see, 
there are important properties of the solutions which do depend on those masses. 

Further properties of the solutions and certain phenomenological classification 
of non-universal boundary conditions follows from the equation (HD [[HJ. Since 
m 2 = m H 2 + /A this equation can be interpreted as an equation for /i 2 . In the 
universal case large values of /i 2 are needed to cancel large negative values of m 2 H2 . 
Now, with non-universal scalar terms it follows from eqs.(§) and (0) that 

/x 2 = c M M 1/2 + c m 10 2 16 + 2D - + . . . , (17) 

where 

9 \ A 3 \ m 2 - m 2 6 



7 yz J V 7 yz J m\ + m 2 6 



and the second term in eq.(|l8) is generated by departures from universality The 



values of /z 2 depend, contrary to equation (16), on the pattern of the deviation from 



universality in and m w . We obtain the following classification |T0 
(A) \i 2 > c M M x , 2 for c m > 



T 2 



(B) /j, 2 < c M M 1/2 for c m < 
It is clear that in case (A), generically the values of fi remain large fi 3> Mz even 
in the limit M 1( / 2 — 0, due to the positive correlation with the scalar masses. 

In case (B) the parameter /i can be arbitrarily small due to the cancellation be- 
twen the scalar and gaugino contributions in eq.(|T7D (with c + 4d > and only then; 
in the opposite case Mi/ 2 > mo, as in eq . (|15|) , and the cancellation is impossible 
unless — c m > cm which is very difficult to achieve). Note that for phenomenologi- 
cal reasons (experimental bounds) we are actually not interested in the strict limit 
Mi/2 = 0. Thus, radiative breaking can be driven by m ^ M 1 / 2 ^ \i — O(Mz). As 
shown in refs.|l0|, UlTj , it is the case (B) which is phenomenologically acceptable, 



with higgsino-like lightest chargino and neutralino. The non-universalities of type 
(A) suffer from similar problems as the universal case and are disfavoured by the 
combinations of constraints from Qh 2 < 1 and BR(b — > s , y). 

The condition c m < puts non-trivial contraints on the values of mi e and fn w 
at the GUT scale and on the Yukawa coupling. In the following we demonstrate 
that they are satisfied by the values of the masses obtained by RG running in the 
minimal SO (10) model from the Planck scale, with universal boundary conditions 
Mi/2 and mo for the soft gaugino and scalar masses at the Planck scale (the unbarred 
quantities denote Planck scale parameters) and for interesting range of values of the 
top quark Yukawa coupling. 

The set of the relevant SUSY-SO(IO) RG equations is as follows: 

^j-a = -ba 2 , (19) 
iy= fHa-UrV, (20) 



dt V 2 
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-A=—aM-UYA, (21) 
dt 2 v ; 

d , V 

—M = -baM, (22) 
^m? 6 = -5Y(m? + 2m? 6 + A 2 ) + ^aM 2 , (23) 
—m 2 10 = -AY(m\ Q + 2m 2 w + A 2 ) + 18aM 2 . (24) 

(JjL 

Here i = ^log^p, « is the SO(10) gauge coupling, M is the running gaugino 
mass, A - the trilinear soft breaking term and as earlier Y = h 2 /4n is the Yukawa 
coupling for the third generation. In eqs.flIDp-(piD we have explicitly introduced 
the numerical values for the /^-functions which depend only on the representation 
assignment of the matter and light Higgs fields. The coefficient b in eqs.(ID) and (22) 



depends on the heavy Higgs field content of the model. We do not have to specify 
it here as our results are stable under varying b in the range from 3 to the maximal 
value of order 30, corresponding to the Landau pole for the gauge coupling^]. The 
solution to those equations read: 

a = TTTT-^ M = M 1/2 -, Y = Y f° , (25) 
1 + ba t a 1 + 14Y F G 



Eg= (?) ' Fg = J E ^ dt 



63 
2 



+ ba 



0,^(63/26+1) ■ 

a ) 



(26) 



The parameters a , M 1 / 2 and Y are the Planck scale boundary values of the gauge 
coupling, gaugino mass and Yukawa coupling, respectively. 

Similarily as for the running from Mgut to Mz it is convenient to define the 
parameter 

Vg = 77 (27) 

YfG 

where 

v r vc\/r \ E G (M GUT ) 

Y fG = hm Y (Mqut) = , r • (28) 

y ->oo (Mgut) 

We can then express the parameter y z introduced in eqs.([10|) and (|1^) in terms of 

VG- 

yz = t- 29 

l + xy G 

where 

x = 7F MSSM (M z ) Y fG . (30) 

The value of x in eq.fl30|) depends on the scales Mgut and Mpi, on the value of the 
gauge /3-function coefficient b and on a^(Mz) (we take here the attitude that Mgut 

■f In models where the heavy Higgs sector correctly breaks SO(10) to the standard model gauge 
group without any additional unwanted massless fields, the beta function is typically b > +7. We 
use b = +3 for numerical analyses in this letter as a conservative choice. The generated non- 
universality is larger for larger values of b. However, as we will describe below, the final result is 
rather insensitive on the value of b. 
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is determined by the crossing of ai and a 2 and we allow for a small mismatch of 
at that scale [0). For M GUT = 2-10 16 GeV, M Pl = 2.4- 10 18 GW, b = 3 and a 3 (M z ) 
in the range 0.11-0.13, we get x in the range 22-25. It increases (decreases) by 5 for 
b = +30 (for Mpi/Mgut larger by factor 10). The condition c m < 0, eq.fllSD, now 
reads: 

2XVG ~ 7 (31) 



mj - m\ % 



"w + m ie 10xy G + 7 
The solutions c m = are shown in Fig.l as the solid curves, for three different values 
of x = 15, 20, 25. Fig.l illustrates the interplay between the GUT scale values of the 
scalar soft masses and the Yukawa coupling which is necessary to assure c m < for 
different values of x. 

As the next step, we solve the eqs.(f2l 



> 



and (|24D and get for the scalar masses: 



10 



12 



m 



m 



16 



15 
14 3 



2/GP0 + ya 1 



where 



and 



a 2 

«3 



- Vg (1 - Vg) 
+ Vg (1 - Vg) 



a x Al + a 2 A M 1/2 + a 3 M? /2 



1 /an\( 63 / 2fe ) 



(32) 
(33) 

(34) 



1_ (-o) {63/2,,+1) & 



63 
26 



v a J 



- 1 + J/G 



/ X (63/26-1) 

63 / 63\ l-(^) 

^26 V 267 t _ /ajA (63/26+1) 




" 26 V 267 1 _ ^aoy DJ 
.267 6 V 26/ ! _ («oy 63/2fe+1) 



(63/26) 



1 _ (ao) (63/26+1) 



(35) 



With explicit solutions eqs. (|32|) and ( p3|) we can check if the relation ( |3lD is 
indeed fulfilled in the model. The dashed lines in Fig.l show the solutions (]32f) 
and ( p3|) for three different values of the ratio Mi/2/mo at the Planck scale and 
with the values of the other relevant parameters as specified above eq.([H]), with 
as(Mz) =0.11 and Aq = 0. It is clear that the solutions to the RG running from 
the Planck scale to the GUT scale satisfy the constraint c m < (solid lines) for 
values of yo close to the quasi-IR fixed point or lower than about 0.2. 

This discussion nicely illustrates the role of the boundary values at Mqjjt for the 
scalar masses in obtaining solutions to radiative breaking in the MSSM with small 
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fi values. However, since those values depend on both Planck scale parameters m 
and Mi/2, it is more convenient to rewrite eq.(17) directly in terms of the Planck 
scale parameters for effective study of the parameter space. Using eqs.(32)-(35) we 
get: 



/< 2 



c M + ^ (3yz - 2) a 3 



M 2 1/2 + \[{2- 3y G ) - 9 (1 - y G ) (1 - y z )\ m\ + ID + 



7 (36) 
where the dots stand for A dependent terms. The M\m coefficient remains always 
positive. The D term must be positive (see eq.(13) and (14)) and in principle, 
with D = drriQ, it should be included into the coefficient. However acceptable 
solutions to radiative breaking are obtained already with very small values of d, 
of order (9(0.01)|. At this point it is worth noting that the negative numerical 
coefficient c in eq.(14) (which is obtained from numerical integration of the 1-loop 
RG equations) goes strictly to zero for y G — > 1. This result follows from the structure 
of the RG equations and explains why very small positive d is sufficient to change 
the pattern of solutions into those of eq. (|16|) (the Aq contribution to eq. (|14]) is small 
but also positive). 



Thus, the necessary and sufficient condition for cancellations in eq. (|36|) to be 
possible is the negative sign of the tuq coefficient: 

> [(2 - 3y G ) - 9 (1 - y G ) (1 - y z )\ . (37) 

By using eq. (^) we easily see that this coefficient is always negative for x in the range 
(0.6-14.4). However, for the values of x generic for the minimal SUSY-SO(IO) we 
obtain non-trivial constraints on the value of the Yukawa coupling. For x = 22 the 
eq.(|37j) is satisfied for y G < 0.2 or y G > 0.6, in agreement with the results presented 
in Fig.l. 

To study in more detail the x dependence of our result (or equivalently, for fixed 
values of M G ut and Mpi, its dependence on a 3 (M^)) and its sensitivity to two-loop 
corrections in the RG running of the gauge and Yukawa couplings below M G jjt, 
we plot in Fig. 2 our two-loop numerical results as a function of a 3 (M z )- Values 
of y G above and below the band depicted by solid lines satisfy eq.(^). For easy 
interpretation we also plot the curves of constant top quark pole masses. One can 
see that for a^iMz) = 0.11 the top quark has to be heavier than 181 GeV or lighter 
than 170 GeV, both regions being of phenomenological interest. For a^iMz) = 0.12 
both bounds move up by about 5 GeV. 

Finally, we comment on the prameter space which gives correct radiative breaking 
and on the sfermion masses. For instance, for m t = 182 GeV, a^iMz) — 0.11 and 
d = 0(0.01) we get /j 2 = D(2)Mf /2 - O(0.05)mg - O(0M)A 2 and e.g. for solutions 
with /i ~ 100 GeV and with fx/Mi ~ 1 we need mo/Mi/2 = 0(5)[J]. In this simple 

§ The maximum possible value of d can be calculated if one specifies the model completely. 
For instance in the models where the rank is reduced by Higgs fields in 16 and 16 representations, 
one calculates the difference in their soft masses m 2 and m 2 using RG equations and obtains 
D = (m 2 — m 2 )/10. m 2 and m 2 can differ easily by 0(1) because of large group theory factors in 
SO(10), and hence d as large as O(0. 1) is possible. 

^ Small values of aziMz) are obtained from the fits to the electroweak data in the MSSM jl6| 

II Contribution from the Aq term can lower mo somewhat. 



7 



example the particle spectrum contains a light pseudoscalar and a higgsino-like 
chargino, both with masses below M z and within the reach of the Tevatron and 
LEP2. It is interesting to note that similar spectrum is predicted from the best fit 
to the electroweak data in the framework of the MSSM |TH . 



The sfermion masses for the third generation tend to remain relatively small, 
too. Combining eqs.(|32|), ([[||) with (||) and (|) we get: 

™q,u,d = ^ [(2 - 3y G ) + 12 (1 - y G ) (1 - y z )\ m\ + c Q ^ D Mf /2 + ... 

~ O{0M)m 2 + 0{A)Ml /2 + . . . (38) 

and analogously for the sleptons, with cl\ E ~ 0(0.2). In our example, the masses 
of the third generation squarks are in the range 200-300 GeV and for sleptons they 
are 0(100 GeV). A more detailed study is necessary to check if their contribution 
to the breaking of the custodial SU(2) symmetry is consistent with the electroweak 
data. 

Summary: In SO(10) SUSY-GUT, a specific pattern of non-universality in the 
scalar masses at the GUT-scale is generated by the RGE evolution from the Planck 
scale to the GUT scale and the D-term contribution induced by SO (10) breaking. 
This particular pattern of non-universality can make the t-b-r Yukawa unification 
phenomenologically viable, consistent with correct electroweak symmetry breaking, 
experimental and cosmological constraints. The top quark mass is either below 
170 GeV or above 181 GeV if a 3 (M z ) = 0.11. 
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FIGURE CAPTIONS 



Fig. 1. Ratio of masess (mf — mf 6 ) / (wf + rn\ 6 ) at the GUT scale as a function of 
the parameter x/q defined in eq . fl27|) . The solid curves A, B and C are the 
solutions to the condition c m = for three different values of the parameter 
x = 15, 20 and 25, respectively. The dashed curves represent results in the 
SO(10) model for three fixed values of the ratio M 1 / 2 /m = 0, 0.5 and 1.0 
(curves a, b and c, respectively) at the GUT scale. 

Fig. 2. The region in the yG~ a ?,(Mz) plane (outside the band between the solid curves) 
in which the condition fl57p can be satisfied in the minimal SUSY-SO(IO) 
model. The deshed curves correspond to the fixed values of the top quark 
pole mass in GeV. The results were obtained by integrating numerically the 
two-loop RG equations below Mqut- 
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